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O . Abstract 
(N 

^ ' We study how present data probe standard and non-standard properties of neutri- 

Q ■ nos and the possible existence of new Ught particles, freely-streaming or interacting, 

among themselves or with neutrinos. Our results include: '^u < 0.40 eV at 99.9% 
C.L.; that extra massless particles have abundance AiV^ = 2 it 1 if freely- streaming 
^ \ and AA*"^ = ± 1.3 if interacting; that 3 interacting neutrinos are disfavored at 

\Q • about 4a. We investigate the robustness of our results by fitting different sub- 

^ ■ sets of data. We developed our own cosmological computational tools, somewhat 

^ . different from the standard ones. 

O 

O 

CLc 1 Introduction 
6 

' Thanks to recent data about the Cosmic Microwave Background (CMB), Large Scale Structures 

■ (LSS) and also Type la Supernovse (SNe), cosmology has become the most sensitive probe of some 
neutrino properties (e.g. of neutrino masses: oscillation experiments test squared-mass differences, 
and other means of probing the absolute neutrino mass are currently less sensitive) and a very 

^ ' sensitive probe of other neutrino properties, including non standard ones [1, 2, 3]. In this paper 

■ we study how present cosmological data determine standard and non-standard 'neutrino cosmology'. 
This includes three different issues. 

i) testing neutrinos: their masses, abundances, . . . 

ii) do photons, neutrinos and gravitons make up the complete list of light particles? Data from 
particle physics allow extra light particles that are neutral under the Standard Model (SM) 
gauge group, and such extra light particles appear in many speculative extensions of the SM, 
one interesting example being simpler string models.^ 

iii) The two above issues can be connected, because neutrinos are the least tested light particles 
and can easily interact with new light neutral particles, in a way that affects the evolution of 
cosmological inhomegeneities. 



^Within the string scenario light particles can be avoided at the price of assuming that strings vibrate on comphcated 
enough higher dimensional geographies [4], such that predictivity seems lost. 
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Table 1: Schematization of the cases considered in this paper. For each one we list the notation of 
the relevant parameters probed by cosmology and refer to the relevant part of the text. 

In section 2 we characterize the fundamental theories and describe the cosmological parameters that 
we want to extract from present data. Since our implementation of the cosmological computational 
tools needed for this analysis somewhat differs from the standard one, we describe it in section 3. 
Section 4 describes our results (table 1 might help in navigating the paper), summarized in the 
conclusions. 

2 Theory 

The effects of new light particles on the evolution of cosmological inhomogeneities are often presented 
in terms of "non standard neutrino properties" because i) particle physics has tested neutrinos less 
thoroughly than other particles, leaving room for surprises; ii) neutrinos z/ are the particles that can 
naturally interact with some new light neutral states. In this section, we discuss how established 
data and theory restrict the behavior of possible new light states, and in the rest of the paper we will 
consider scenarios that are compatible with these restrictions (more exotic possibilities are sometimes 
entertained in purely phenomenological analyses). 

We consider light particles, neutral under all SM gauge interactions. Indeed the LEP measurement 
of the invisible Z width [5] implies that hypothetical new light particles must be neutral under SU(2)l 
gauge interactions and can have at most a small hypercharge; furthermore new light particles with 
strong interactions seem excluded.^ Field theory and gauge invariance significantly restrict the 
interactions of new neutral light particles in a well-known way that depends on their spin. 

• New light fermions, neutral under all SM gauge interactions (commonly called 'sterile neutrinos' 
Vs), can have a mass mixing with ordinary 'active' neutrinos. This is described by the following 
SU(2)l U(1)y gauge-invariant Lagrangian, written in compact notation 

^{u,) = ^sM + i^si^i^s + {XLHus + y + h.c). (1) 
L and H are the lepton and Higgs weak doublets. 

• New light neutral scalars (p can have SU(2)l ® U(l)Y-iiivariant Yukawa interactions with sterile 
neutrinos, as described by the following Lagrangian: 

^(z^s, cP) = ^{v,) + {d^cl^f + (^ <^ + h.c.) + y(0). (2) 

^On the experimental side, data can leave open windows. On the theoretical side, light colored particles typically 
form hadrons not lighter than the QCD scale; one can however imagine a colored scalar with a negative 'tachionic' bare 
squared mass, fine tuned to almost cancel the QCD energy. We name this speculative particle 'rattazzon' in honor of 
the only theorist that accepted to discuss it. 
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Table 2: Differences between our approach and a 'typical' cosmological analysis. 

At energies below the mass of a sterile neutrino, it can be integrated out so that one effectively 
obtains Yukawa couplings vi'(j), suppressed by the mass of the sterile neutrino, as well as Yukawa 
couplings i^fs0 with other lighter sterile neutrinos. We are aware of no reasonable allowed way 
of coupling a scalar directly to active neutrinos, in absence of light sterile neutrinos. 

• Sizable couplings of neutrinos to light vectors can be introduced in a similar way: adding sterile 
neutrino(s) charged under an extra spontaneously broken gauge symmetry. 

The gauge and Lorentz symmetries allowed theorists 30 years ago to successfully predict and guess 
experimental results. Theorists have tried to proceed further by demanding a further restriction: 
naturalness. At the moment it is unclear whether this is a correct requirement; for instance, the 
LHC will tell whether physics at the weak scale obeys it or not. Its imposition would strongly 
restrict the behavior of new light scalars. For example, one could explain their lightness by assuming 
that they are Goldstone bosons of spontaneously broken lepton numbers: in models with a single 
scalar, this implies that it couples to neutrino mass eigenstates (rather than to generic combinations), 
forbidding neutrino decay in vacuum. 

Even taking into account that these considerations significantly restrict the set of 'reasonable' 
fundamental theories, flavour issues make the number of fundamental parameters so large that a 
direct study of the fundamental theory seems unpractical. (This approach has been pursued in [6] in 
the case of a single sterile neutrino and no bosons). On the other hand, cosmology probes (new) light 
particles via their gravitational couplings: neutrinos and new light particles affect the evolution of 
inhomogeneities in the gravitational potential(s), in which matter moves. Therefore it is convenient 
to focus on parameters half-way between theory and data: equation of state, sound speed, etc. We 
will consider the following basic limiting cases: extra freely-streaming particles, or extra particles 
that interact among themselves, or extra particles that interact with neutrinos. Whether a particular 
model falls in one or the other category depends on the value of the fundamental couplings (like A, 
A' and the interactions in V{(j)) above) of the particular model. 

In the following analyses we will study different systems specifying the parameters probed by 
cosmology. Some of these parameters will be combinations of the parameters appearing in the low- 
energy effective Lagrangian, some others will not (e.g. the abundances around recombination). The 
latters could be translated in terms of the fundamental parameters of each full model, a step that 
we do not make. 

We ignore other possible probes of the same physics (which e.g. also manifests as neutrino decay, 
anomalous matter effects, rare K ^ iv decays), as they presently are orders of magnitude less 
sensitive than cosmology. 
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Figure 1: Our computation of CMB and LSS spectra in standard ACDM cosmology, compared with 
data. 



3 Analysis strategy 

We develop our own computational tools for the analysis of the cosmological observables. For what 
concerns standard cosmology, our results agree with those of other authors (e.g. the WMAP Sci- 
ence Team), but having independent analyses is clearly important. In this respect, our analysis is 
particularly independent: it differs from what nowadays is a typical analysis in the way illustrated 
in table 2. Cosmological observables are computed using a code written by one of us, rather than 
running the commonly used CMBfast or CAMB public codes [8]: this allows us to have a better 
control and flexibility on non-standard modifications. 

We use the line-of-sight approach in the conformal Newtonian gauge [8, 9, 10]; recombination 
can be implemented both in Peebles approximation (see e.g. [31]) and using the external recfast 
code [11], which is the option chosen for the present analysis. 

The main disadvantage is that our code is almost 2 orders of magnitude slower than CMBfast 
or CAMB. In part this happens because, rather than optimizing our code for standard cosmology, 
we keep it fully flexible such that non-standard cosmologies are immediately implemented.^ In 
part this happens because, while standard codes are written in FORTRAN, our code is written in 
Mathematica [12] and we run it on a common laptop (rather than on a cluster of computers). 

We now describe the advantages of our approach that allowed us to perform our analysis. Readers 
not interested in these technical details can skip the rest of this section. The main point is that, 

^For example, the interacting particles considered in this paper are implemented by typing their linear evolution 
equations, eq. (10) or eq. (13), in NDSolve form. 
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Figure 2: Difference between our code and CAMB, at the standard- cosmology best-fit point for m^j = 
(red solid line) and form^ = 0.5 eV (blue dashed line). Our code does not employ any approximation 
specific for standard cosmology. In both codes various parameters allow the user to increase the 
accuracy; this plot holds for the choice employed in the present paper. The dotted line shows the la 
accuracy obtained by WMAP3 results (binned data), indicating that we have a good enough accuracy 
(as confirmed by other tests). A similar %-level accuracy is found for the TE and EE CMB spectra, 
that presently are measured with much larger uncertainties than the TT spectrum. 

while FORTRAN can only do numerical computations, Mathematica does not have this limitation and 
allows to do analytically all parts of the computations that can be done analytically. This includes 
the dependence of cosmological observables, e.g. on the spectral index, and all statistical issues that 
nowadays are the most time-consuming aspect of cosmological analyses. Our approach is based on 
the powerful old-fashioned Gaussian techniques, as we now briefly describe. 

3.1 Statistics 

Cosmological data have become so accurate and rich that debates about Bayesian priors versus fre- 
quentistic constructions are getting numerically irrelevant: all different techniques converge towards 
their common gaussian limit. This is clear e.g. from figures 10 of the WMAP analysis [7]: within 
good approximation all allowed regions identified by their Monte Carlo Markov Chain (MCMC) tech- 
nique are ellipses (with sizes that have the Gaussian dependence on the confidence level), as they 
must be in Gaussian approximation. This means that the usual x^, a single quadratic function of 
the various cosmological parameters, approximatively encodes all present information on standard 
cosmology and that the dependence on the iV^*"'^ parameters of standard cosmology (here chosen 
to be the usual As,ns,^-DM-,^h-,yp-,h,T with Otot = defined as in [7]; Ag is normalized at the 
pivot point k = 0.002/Mpc) is accurately enough described by a first order Taylor expansion of each 
observable (the power spectra, the luminosity distances of supernova, 

...) around any point close enough to the best-fit point. We will soon check explicitly that sampling 
Nf""^ + 1^10 points is enough to study standard cosmology.'^ For comparison, MCMC techniques 

*We do not improve the accuracy by making a second-order Taylor expansion. This can be done by probing 
jYstnd|-^stnd ^ ]^y2 < 50 morc points only, but would complicate statistical issues, preventing e.g. analytical marginal- 
izations of the likelihood over nuisance parameters. 

Furthermore, we checked that using two-sided derivatives or recomputing observables with the public CAMB code [8] 
affects the results of the global standard fit, eq. (3), in a minor way: best-fit values typically shift by about a third of 
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Figure 3: Fit of cosmological data at 68, 90 and 99% C.L. The shaded areas show our global fit 
without Lyman-a, and the dotted lines our WMAP3-only fit, such that this figure can be directly 
compared with the analogous WMAP Science Team plots in fig. 10 of [7j. 

can need one independent chain with ~ 10^ points every time one wants to analyze a (sub)set of 
data [7]. 

The Gaussian approximation has no 'statistical' uncertainty due to finite MCMC sampling but it 
introduces a 'systematic' uncertainty. This is small near the expansion point (chosen to be close to 
the best-fit point) and grows when one goes far from it. At some point data become accurate enough 
that the region singled out by them is small enough to make the Gaussian approximation a good one. 
By construction, the Gaussian approximation reproduces the same best-fit point (small differences 
between our and other analyses on common studies are due to different data-sets, different code, etc.) 
and the confidence regions with small enough confidence levels, and fails at larger confidence levels. 
In practice, we care about 90%, 99% and maybe 99.9% confidence levels. Fig. 3 is our crucial test 
and it shows that the contours corresponding to such confidence levels are reproduced in an fairly 
accurate way when comparing with the WMAP Science Team analysis. Notice that the Gaussian 
approximation needs not to be and is not accurate enough to analyze every single piece of data, 
but it allows to correctly fit the full data-set. Some non-standard cosmological parameters are still 
subject to 'degeneracies': we will later discuss how the Gaussian approximation can be extended to 
deal with these situations. 

Our code directly gives the ^is an analytic quadratic function of cosmological parameters, that 
fully describes present information on ACDM cosmology. Our result in terms of best fit points and 
la errors is 

fit As h Us T lOOf^fe/i^ ^TiMh"^ 

WMAP3 0.80 ±0.05 0.704 ± 0.033 0.935 ± 0.019 0.081 ± 0.030 2.24 ±0.10 0.113 ±0.010 (3) 
Global 0.84 ±0.04 0.729 ± 0.013 0.951 ± 0.012 0.121 ± 0.025 2.36 ±0.07 0.117 ±0.003 



a standard deviation. 



The symmetric correlation matrix for the global fit is^ 

/I - - - - - \ 

-0.19 1 

-0.25 0.71 

^ ~ 0.74 0.31 

-0.08 0.70 

0.07 -0.03 0.14 -0.02 0.36 1 / 

Fig. 1 shows our best-fit, and fig. 2 shows the percentage difference in the main observables among 
our code and CAMB [8]. 

3.2 The data set 

Our global fit includes the following set of data: 

• WMAP: the 3 year TT, TE and EE data from WMAP [13] (called WMAP3 throughout). 
We use data in binned form (72 bins) and deduce the Gaussian approximation of the full 
likelihood by extracting mean values, errors and their correlations from the numerical likelihood 
code provided by the WMAP collaboration [14], that includes the uncertainty due to cosmic 
variance. At small t there is a non-Gaussian large uncertainty; a fortiori these data do not 
have a dominant statistical weight. We verified that using data in unbinned form makes a very 
minor difference. 

• Other CMB data (103 bins): the most recent data from ACBAR, Boomerang (TT, TE and 
EE results), CAPMAP (EE), CBI (TT and EE), DASI (TE and EE), VSA [15]. 

• Large Scale Structures: 

— The SDSS [16] and 2dF [17] measurements of the matter power spectrum from galaxy 
surveys (called LSS throughout). An important and long-standing issue in this respect 
is how to relate the quantity of interest for cosmology (the matter power spectrum P{k), 
computed with linear perturbation theory) to the measured galaxy-galaxy power spectrum 
-fgai(fe)- We adopt the prescription suggested in [17] to effectively take into account the 
normalization bias and the nonlinear effects: -Pgai(^) = ft^(l + Qfc^)/(1 + Ak)P{k). We 
adopt A = lA and we marginalize over a free bias parameter b and over Q = 4.6 it 1.5 and 
10±5 for 2dF and SDSS respectively (we assume Gaussian errors), following [17] and [18]. 

— The detection of Baryon Acoustic Oscillation peaks (BAO) in the correlation function of 
the SDSS subsample of Luminous Red Galaxies at nominal redshift z = 0.35 [19]. We 
implement it in terms of a measurement of the adimensional parameter 

A = D,{0.35)^n^^tterHl/{0.35 c) = 0.469 ± 0.017 

This is closely related to the Fisher matrix widely used by cosmologists. Since diflterent communities use different 
terminologies, we here define the precise meaning of the quantities we employ. The mean values fii of the parameters 
Pi, their errors ai and the correlation matrix ptj determine the a-s 

= ^{Pi - fJ■^)[('^^^^h(PJ - W), where {a'^),j = cnpijaj . (4) 

id 

The likelihood is given by £ = £maxe^^ . In Gaussian approximation, marginalizing £ with respect to any sub-set of 
'nuisance' parameters is equivalent to minimizing the with respect to nuisance parameters. E.g., since eq. (4) contains 
the inverse of the matrix, the marginalized over all parameters p except one of them, pi, is X^{Pi) = {Pi ~ I^O^/'^i- 
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(neglecting a small dependance on the primordial spectral index Ug), where c is the speed 
of light and is a distance defined in terms of the comoving angular diameter distance 
Da and the Hubble parameter H as D^{z) = \Dx{z)'^cz/ H{z)\^^'^ . We computed the 
galaxy-galaxy correlation function in a few cases, checking that this is a satisfactory ap- 
proximation. Within standard cosmology, BAO data give a measurement of the total 
matter density: 

^matter = 0.28 ±0.025 [19]. (6) 

— The Lyman-a Forests in distant quasars absorption spectra (called Lyman-a through- 
out) from Croft et al. [20], at redshift z = 2.72, and from SDSS as condensed by [21] 
into the measurement of the renormalized amplitude of the power spectrum /\^{k,z) = 
k^P{k,z)/{2TT'^) = 0.452 ± 0.072 and of its slope nes{k,z) = din P{k) /dink = -2.321 ± 
0.069 at the pivot point k = 0.9 /i/Mpc and z = 3.0.® Different re-analyses of the same 
data find somewhat different results (see e.g. [22, 18]), which might signal systematical 
problems, e.g. in the way the flux power spectra are converted into measurements of the 
matter power spectra. In the following we will adhere to the values presented in the data 
sets above, but pay special attention to the implications of their use. 

• Type I supernovas: the Gold sample of Riess et al. [23] and the SNLS05 data [24]. We combine 
the datasets directly, although strictly speaking they are not independent as they share the same 
set of low-redshift supernovse and use slightly different techniques in the analysis; however these 
simplifications marginally affect the final results such that a fully careful combination seems 
unnecessary; present cosmological data surely contain more worrisome issues. Within standard 
cosmology, supernova data give measurements of the total matter density, 

^^mattcr = 0.33 ± 0.036 [23], n^^tter = 0.26 ± 0.038 [24]. (7) 

• The HST measurement of the Hubble constant Hq = 100 /ikm/sec/Mpc: h = 0.72 ± 0.08 [25]. 

All CMB data are marginalized with respect to the recombination and Sunyaev-Zeldovich back- 
grounds. 

Our strategy concerning statistics needs to be extended to study the new physics scenarios ad- 
dressed in the next section. In general the dependence of the on the new physics parameters is not 
Gaussian, because we consider small new physics effects that behave in significantly different ways 
in different regions of their parameter space. 

In most cases (e.g. massive neutrinos, section 4.1) the simplest extension is enough: we compute 
the observables as function of the new physics parameters (e.g. for a dozen of values of neutrino 
masses) and, at each point, we use the Gaussian technique to marginalize over the standard cosmo- 
logical parameters listed in eq. (3). 

In other cases (e.g. in presence of extra massless neutrinos, section 4.2) this turns out to be 
not enough, because data still allow for large new physics effects: e.g. a large variation in Nu can 
be compensated, to a large extent, by readjusting the other cosmological parameters, expecially h 
and r^DM [29]. We deal with this kind of 'degeneracy' by applying the well-known Newton iterative 
minimization technique: for any given value of the new parameters (e.g. N^) we perform the Gaussian 
expansion around some point: in general it gives an inaccurate value of the minimum because 

®We reduced the errors by 25% because this brings the Gaussian in closer agreement with the 'exact' hkehhood 
as provided in [21]. 
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Figure 4: Fit as function of the energy density in freely-streaming slow particles, parameterized by 
the sum of active neutrino masses, assumed to have the standard abundance = 3.04. We studied 
different combinations of data-sets, as indicated by the legend. 

the expansion point is not close enough to the best-fit point: we then use the just estimated best-fit 
point to perform a second Gaussian expansion around it, finding an improved approximation to the 
minimum and to the true best-fit point. In practice, two or three iterations of this procedure are 
enough to find a good approximation, as can be confirmed by checking that more iterations no longer 
change the result. Subsets of data can then be analyzed by performing the Gaussian approximation 
around this global best-fit point. 

4 Results 

In all cases we will report the difference Ax^ with respect to the standard ACDM model (which is 
individuated by the best fit points in eq. (3) and assumes the existence of only 3 ordinary, massless 
neutrinos), because this is the significant statistical indicator. Roughly speaking, Ax^ + ANp = itn^ 
corresponds to n standard deviations of evidence against or for, where AA'p = Np — Np^^'^ is the 
number of extra relevant free parameters with respect to ACDM. For the models that we will consider, 
the not univocally defined AA'p ~ 1 or 2 is small enough that it does not play a significant role. In 
practice, this means e.g. that 3a ranges can be read from fig. 4 or fig. 5 by looking where our curves 
reach Ax^ = 9. 

We here do not judge the compatibility of a model with data by performing the Pearson test. 
According to this test, a model is excluded if the total x^ is not compatible with its expected value 
A^ lb \/2N where A^ is the number of data points minus the number of free parameters (for A^ ^ 1 the 
X^ follows a Gaussian distribution with this mean and variance; for small N the x^ has a different 
'x^' probability distribution). The \/2N means that the Pearson x^ t^st becomes inefficient when 
N ^ 1, which is the case in cosmological fits, where A^ ~ (100 1000). In practice, this means that 
a model disfavored at 5a by the Ax^ test typically is disfavored only at la by the x^ test. A more 
detailed discussion of these issues can be found in section 4 of [26] (in the context of fits of solar 
neutrino data). 
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4.1 Massive neutrinos 



Fig. 4 shows the resuh of our comparison of cosmological data with massive neutrinos. We here 
assume that only the three active neutrinos exist, and have the standard abundance and temperature. 
Cosmology is not yet sensitive enough to discriminate neutrinos with normal from inverted hierarchy, 
but only to heavy quasi-degenerate neutrinos. Therefore the key parameter probed by cosmology 
is the sum of neutrino masses J2''^u = fni + + fn?,- Oscillation data imply '^rrii, ~ 0.05 eV 
(0.10 eV) [1] if neutrinos have strictly normal (inverted) mass hierarchy, and larger values of "^iriy 
are obtained for quasi-degenerate neutrinos. We assumed normal hierarchy at J2^i> — 0.05 eV, 
inverted hierarchy at "^i^ = 0.1 eV, and degenerate neutrinos at larger ^ttt-i,. 

Within standard cosmology, WMAP3 data alone imply a constraint not plagued by potential 
systematical problems: we find ^2^^ < 2.2 eV at 95% C.L. (see fig. 4) in agreement with [7, 27]. 
Adding LSS data (that are more strongly affected by neutrino masses than CMB data) gives stronger 
constraints, but they are subject to the potential systematic problems discussed in section 3.2. We 
therefore plotted different lines, that correspond to different combinations of data-sets, with each 
kind of observable excluded in turn. Our global fit gives 



in good agreement with [28] and in reasonably good agreement with [18] . We see that SN data have 
essentially no effect, while Lyman-a data play an important role: the constraint becomes about 2 
times weaker if Lyman-a data are dropped: "^miy < 0.73 (0.52) eV at 99.9% (95%) C.L., in good 
agreement with [28] and in reasonably good agreement with [7]. 

However, the constraint from the global fit is somewhat stronger than the sensitivity of the data. 
Indeed, focussing on Lyman-a data, the main effect of neutrino masses at fixed values of other 
cosmological parameters, consists in reducing the adimensional A^, measured by SDSS Lyman-a 
data with ±0.07 uncertainty, by about O.lS^m^/ eV. However, for massless neutrinos, the global 
fit already suggests a value of about 2 standard deviations below the experimental value, and 
making neutrinos massive reduces further the expected A^. In eq. (8) we avoided reporting a strong 
but doubtful constraint by choosing a confidence level higher enough than 2 standard deviations. 

4.2 Extra freely-streaming massless particles 

As usual the 'number of neutrinos' N^, parameterizes the amount of energy in all relativistic freely- 
streaming degrees of freedom, converted in terms of 'neutrino equivalents': Ni, includes the ordinary 
neutrinos and any extra fermion or boson, as it is formally defined by the relation /Orelativistic = 
[1 + 7/8 iV,.(r^/r)4] , where T^/T = (4/11)^/3 at T < me. Standard cosmology with 3 neutrinos 
predicts Ni, ~ 3.04 (the deviation from 3 being due to the incomplete decoupling at the time of 
the beginning of e^e~ annihilations, plus other small corrections). Our global fit gives 



and the continuous line in fig. 5a shows the precise form of the A^^. We checked that using the first 
release of WMAP data our result agrees with the various analyses published in the literature [30] . 

In general, results of global fits can be mislead by problems in any of the pieces of data they 
contain; in this case the validity of the global fit appears particularly doubtful: N^, is dominantly 
determined by non-CMB data, and giving slightly different weight to them can significantly affect the 
fit because different pieces of data prefer different values of N^. In particular, the 2a preference for 



^m^<0.40eV at 99.9% C.L. 



(8) 



N^ = 5±l 



(9) 
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Figure 5: Fig. 5a): fit as function of the energy density in freely- streaming relativistic particles, 
parametrized by the usual 'number of neutrinos' N^. Fig. 5b): fit as function of the energy density 
in extra interacting relativistic particles, with abundance parametrized by AN^. We studied different 
combinations of data-sets, as indicated by the legend. 

Ny > 3 is mainly due to the 2a anomaly in the Lyman-o measurement of the power spectrum: fig. 5a 
shows that omitting Lyman-a one recovers excellent agreement with the standard value Ny = 3. The 
agreement with and between the up-to-date analyses performed by the WMAP Team [7] and by [18] 
is imperfect; in particular the revised version of [18] claims a 3cr preference for Ny > 3. 

4.3 Extra massless particles interacting among themselves 

In the previous section we considered extra (massless) particles with negligible interactions, that 
therefore freely move on cosmological scales. We now consider the opposite limit: extra (massless) 
particles that interact among themselves with a mean free path smaller than relevant cosmological 
scales, such that inhomogeneities in their energy density evolve in a different way. Concrete examples 
are an elementary scalar with a quartic self-interaction or any particle with low compositeness scale, 
obtained e.g. if some extra QCD-like gauge group becomes strongly coupled at an energy much lower 
than the QCD scale. In the tight coupling limit this system is described by a fluid: its density and 
velocity perturbations 5 and v obey the standard fluid equations (in the conformal Newtonian gauge 
and in linear approximation): 

4 k5 

§ = -4$ kv, V = -\ (10) 

3 4 ^ ^ 

where a dot denotes derivative with respect to conformal time, k is the wavenumber, $ and ^ are 
the scalar perturbations in the metric (in the notations of [31]). 

Fig. 5b shows the constraint on the density of the extra particles, that we parameterize in terms 
of the usual 'equivalent number of neutrinos' AA''^^ > 0: the global fit gives 

ANy = ± 1.3. (11) 
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Figure 6: Fig. 6a: how the global fit constrains the presence of extra freely- streaming particles with 
abundance AN^, and mass mg. Shaded regions are allowed. Fig. 6b: the analogous result for extra 
interacting particles. 

Data do not favor the presence of extra massless particles interacting among themselves. Dropping 
Lyman-a data makes the la constraint two times more stringent. Fig. 5b shows that the interval 
allowed at n-a is not n times larger than the la interval. 

4.4 Extra freely-streaming massive particles 

We here assume that the ordinary neutrinos have cosmologically negligible masses (e.g. 0.05 eV [1]), 
and add extra freely-streaming particles, with abundance ANi,, mass nig and the same temperature as 
the ordinary neutrinos. We have chosen the notation ANi, and nis because, although the constraints 
we obtain here apply to a generic freely-streaming fluid (but not to extra interacting particles, studied 
in the next sections), sterile neutrinos of mass rus are by far the most popular specific realization. 
Sterile neutrino models allow one to compute the abundance AN^, in terms of oscillation parameters 
(non-minimal scenarios introduce extra parameters, such as particle/anti-particle asymmetries), and 
ANi^ can acquire a small, sub-thermal value, that stays constant after neutrino decoupling at T ~ 
MeV. We here assume that AA^^, keeps an arbitrary constant value at the much lower temperatures 
relevant for CMB and LSS observables. The present energy density is 

Si. . ^ = «^!^. (12) 

Per h^ eV 

This case is characterized by 2 new parameters: therefore we must make 2-dimensional plots, so 
we prefer not to show how the fit changes by considering various different data-sets. Fig. 6a shows 
the global fit. We can distinguish three regions. For small ttt-s and large AA^^> 1, data disfavor a 
too large number of relativistic particles around recombination. In the intermediate region the extra 
particles behave as warm Dark Matter: their abundance is constrained to be ^ 0.01 dominantly by 
Lyman-a and LSS data. For larger rrig ^ keV the extra particle behave as cold enough Dark Matter 
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and sterile densities as large as Jig ~ 0.22 are allowed. Assuming one thermalized freely-streaming 
sterile neutrino {AN^j = 1) we find that its mass is constrained to be mg < 0.7 eV at 99.9% C.L., a 
range incompatible with the mass suggested by the LSND anomaly [33, 1]. Our results have some 
overlap and substantial agreement with [6, 34, 18, 35]. 

4.5 Extra massive particles interacting among themselves 

We now assume that the extra particles discussed in section 4.3 have a non negligible mass mg 
and are stable, such that when the temperature T falls below nis they form a non-relativistic relic. 
(Alternatively, they could decay into freely-streaming particles realizing a more complicated situation 
that is intermediate between the one studied in the previous section and in this section). Particles 
with this behavior are exemplified in section 4.3. In the tight coupling limit this system is described 
by a fluid, and the massless fluid equations in eq. (10) generalize to 

2 

§ = -(l + yj)(3<l> + kv)-3-(c!-w)5, V = - -n - 3c^)v + -^kS (13) 

a a 1 + w 

where w = p/p and Cg = 6p/6p is the squared sound speed. We fix the equation of state of the fluid 
by assuming that its average energy density p and pressure p is the one of 2 • ANi, dof decoupled 
from the rest of the thermal plasma: 

^ = '-^^^/(0^^^"'''^' ^ = '-^^^/(03|^^-^/^% ^ = VJ^- (14) 

In practice this just means that the fluid interpolates between relativistic and non relativistic matter: 
z/; ~ Cg ~ 1/3 at r ^ nig and ~ Cg ~ at T ^ ttt-s. A more sophisticated treatment seems 
unnecessary. For simplicity we have adopted Boltzmann statistics in eq. (14). Again, the parameter 
AA^,^ tells the initial abundance of the extra particles in the usual 'neutrino-equivalent' units. We 
again assume that Tg = T,^, such that this extra component is described by two parameters: its 
abundance AA*";^ and its mass Trig. 

Fig. 6b shows the result of the global fit: interacting extra particles are constrained in a slightly 
stronger way than freely-streaming extra particles (at not too low C.L.). Notice that, like in the case 
of massive freely-streaming sterile neutrinos, cosmology disfavors the mass values suggested by the 
LSND anomaly [33, 1] also in this opposite limit of tightly-interacting sterile neutrinos, suggesting 
that intermediate cases might also be not viable. 

A. de Gouvea suggested us one possible new economical interpretation of the LSND anomaly, 
in terms of decays among active-only neutrinos. This would need some order one couplings with a 
light scalar, making neutrinos cosmologically interacting: as we have seen a global fit of cosmological 
data disfavors this possibility. Furthermore, a preliminary analysis indicates that the needed decay 
is incompatible with SuperKamiokande atmospheric neutrino data at about 3a C.L. 

^Recent dedicated analyses found that the case in which sterile neutrinos provide all the observed Dark Matter, 
JIdm ~ 0.22, and (as we assume) have the same temperature of active neutrinos, is allowed by cosmological data only 
for ma^lOkeV [32]. (Such large masses are incompatible with X-ray observations, so that this scenario is allowed 
only if these particles are produced with sub-thermal velocities) . We limit our plots to the region nis < 1 keV, where 
the allowed fls is small enough that our Gaussian technique (implemented in the most nai've way) provides a good 
approximation. The structure present in our plot at nis ^ keV might be artificially too sharp, because Lyman-a SDSS 
data have been condensed in measurements of the power spectrum at a single wave-number k. 
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Figure 7: Global fit of the total relativistic energy density (parametrized by the usual 'number of 
neutrinos' N^), and of its interacting fraction. The dot shows the SM prediction, and along the 
dotted line the number of normal freely-streaming neutrinos keeps its standard value. 

4.6 Massless neutrinos interacting with a massless boson 

So far we added extra light particles, free or interacting among themselves. We now assume that 
ordinary neutrinos are involved in interactions with these extra particles. More specifically, we 
consider ]V^°''™'^^ neutrinos that behave normally, while neutrinos are involved in the interactions 
with extra scalars (p, such that these interacting neutrinos no longer free stream, but form a 
tightly coupled fluid together with the scalars. 

Following [36, 37] we assume that the energy and pressure density of this fluid are given in the 
homogenous limit by 

p = 2Nl^'p'^{m„Tfi) + N^p'^im^, Tfl), p = 2Nl^Y^{m,, Tfl) + N^p'^{m^,Tii) (15) 

where p'p^m, T) and p^p{m, T) are the energy and pressure density of one fermionic degree of freedom 
with mass m in thermal equilibrium at temperature T, and and are the analogous quantities 
for one bosonic degree of freedom. The temperature Tg of the fluid is computed by assuming that it 
equals the ordinary neutrino temperature at T ^ rui,, m^, and that the fluid cools down adiabatically. 
This fixes the fluid equation of state w and its sound speed Cg , and inhomogeneities evolve as dictated 
by eq. (13). Summarizing, this system is described by the following parameters: iV"™™'^^ 

In this section we assume that and are negligibly small, such that w = Cg = 1/3 (relativistic 
fluid). Then, the ratio N^^/N^ becomes essentially irrelevant, such that the system can be described 
by just two parameters: i) the total energy density in relativistic particles, that we describe by the 
usual 'number of neutrinos' N^, = J\[^°^^^^ -\- A^™* + 4iV(^/7, that remains constant; ii) the energy 
fraction R = N™^ /Ny that contributes to the fluid. The remaining fraction 1 — R = N^°™'^^ /Ny 
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Figure 8: Number of standard deviations, defined as (x^ — Xo)^^^ (where Xo is the best ACDM fit 
with massless neutrinos) at which cosmological data disfavor a fluid of 3 neutrinos interacting with 
a scalar assuming massive neutrinos and massless scalar (fig. 8a) or massless neutrinos and massive 
scalar (fig. 8b). The different lines correspond to different data-sets: global fit (continuous black line), 
Lyman-a data dropped (dot-dashed green line), WMAP3 only (dashed red line). 



freely streams. In standard cosmology R = and Ni, = jV^™™ = 3.04. Fig. 7 shows how a global 
fit of present data determines these two parameters. The 'all interacting' case {R = 1) is disfavored 
at 4(T at least (i.e. minx^iNj^, R = 1) — x^{Ny = 3,R = 0) > 16) and at 3a if Lyman-a data are 
dropped. As in the case of massive neutrinos, Lyman-a data make the constraint slightly stronger 
than the sensitivity. Two previous analyses claimed different results: our constraints are somewhat 
stronger than in [37] (possibly because we use the most recent data set) and weaker than in [36]. 

4.7 Massive neutrinos interacting with a massless boson 

We now explore how the situation changes if neutrinos have a non vanishing mass m^. We focus 
on the most interesting limiting case: R = i.e. we now assume that all neutrinos are involved 
in the interaction. This is interesting because it means that the cosmological bound on neutrino 
masses no longer applies, because when T < ruy all neutrinos annihilate or decay into massless (f) 
particles. Scenarios of this kind have been proposed for a number of reasons [38, 39, 40]. We again 
assume that neutrinos initially have the standard abundance, and that bosons initially have the 
minimal abundance, N(f, = 1 (one real scalar). After that all neutrinos annihilate into 0, they acquire 
a relativistic energy density corresponding to an equivalent number of neutrinos N^{T <m^) = 
4/7(25/4)4/3 ~ 6.6. 

Fig. 8a shows how much this non-standard cosmology is disfavored as a function of m^, (standard 
cosmology is not recovered for any value of m^). For nii, <^ eV the result is similar to the case 
nil, = 0, already discussed in section 4.6: this scenario is disfavored at about 4cr by the global fit. As 
already noticed in [37], the scenario becomes less disfavored for rui, > eV (beta decay data demand 
777-1/ ^2eV [1]). We find that WMAP3 data (dashed lines in fig.s 8) are more constraining than the 
WMAPl data analyzed in [37]. 

We do not consider intermediate scenarios where only one or two massive neutrinos interact with 
the scalar: both the constraint on neutrino masses and on their free-streaming applies, but in a 
milder form [37]. 
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4.8 Massless neutrinos interacting with a massive boson 

We conclude studying the opposite limit: neutrinos have a negligibly small mass, while the scalar 
has a mass m^j,. The considerations of the previous section still apply, with the role of neutrinos and 
interacting particles interchanged: since neutrinos have more degrees of freedom than one scalar, the 
radiation density increases in a mild way when T drops below m^. Fig. 8b shows how much this 
non-standard cosmology is disfavored as a function of m^: we find almost no dependence on m,^: 
this scenario is disfavored at about 4cr. For large enough m^, depending on the model, interactions 
mediated by (p must become weak enough that neutrinos recover their standard freely-streaming 
behavior. 

5 Conclusions 

We compared a non exhaustive but representative casistics of how cosmology is affected by extra 
light particles (with sub-keV masses), or by standard and non-standard properties of neutrinos, 
using CMB, LSS, Lyman-a, BAO, SN data. 

• First, we considered ordinary massive neutrinos. We obtain the cosmological bound on neutrino 
masses, J2 ''^u ^ O-^O eV at 99.9% C.L. and fig. 4a shows that the relatively less safe observations 
play a crucial role. 

• The density of initially relativistic particles can be parameterized in terms of the usual number 
A^jy of equivalent neutrinos. Assuming that all the Niy relativistic particles freely stream, we 
find that their density is constrained to be = 5 it 1. The 2a preference for Niy > 3 is mainly 
due to the 2a anomaly in the Lyman-a measurement of the matter power spectrum. 

• Assuming ordinary neutrinos plus an extra component of interacting particles, we find ANi, = 
lb 1.3. Fig. 7 shows how data constrain the intermediate case where both kinds of relativistic 
particles are present. It is interesting that the uncertainty on AN^, is decreasing below 1. 

• The extra light particles might have a mass m and an abundance AN^. Fig. 6 shows how data 
constrain these parameters in the two limiting cases that these extra particles freely stream 
(fig. 6a) or interact among themselves (fig. 6b). 

• Finally, we considered one extra scalar of mass m,^ that interacts with neutrinos of mass mjy. 
We find that this scenario is strongly disfavored by the global fit, at about 4(7. 

All these results are based on assumptions and subject to caveats, that we discussed in the text. 
Technically, our analysis somewhat differs from typical analyses because we used a code developed 
by us and dealt with statistics using Gaussian analytical techniques, that become adequate nowadays 
that observations are rich and precise enough. Eq.s (3) and (5) allow to check how well we reproduce 
the standard results for standard cosmology. 
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